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ABSTRACT 

Discovery of the first planetary system by direct imaging around HR8799 has made 
the age determination of the host star a very important task. This determination is the 
key to derive accurate masses of the planets and to study the dynamical stability of the 
system. The age of this star has been estimated using different procedures. In this work 
we show that some of these procedures have problems and large uncertainties, and the 
real age of this star is still unknown, needing more observational constraints. Therefore, 
we have developed a comprehensive modeling of HR 8799, and taking advantage of its 7 
Doradus-type pulsations, we have estimated the age of the star using asteroseismology. 
The accuracy in the age determination depends on the rotation velocity of the star, and 
therefore an accurate value of the inclination angle is required to solve the problem. 
Nevertheless, we find that the age estimate for this star previously published in the 
literature ([30,160] Myr) is unlikely, and a more accurate value might be closer to the 
Gyr. This determination has deep implications on the value of the mass of the objects 
orbiting HR8799. An age around «1 Gyr implies that these objects are brown dwarfs. 
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1 INTRODUCTION 

The discovery of the first planetary s ystem by direct imaging 
around HR8799 |m arois et al . 2008) was an important mile- 
stone in the field of exoplanet research. Up to now, eleven 
planets have been discovered with this procedure, and only 
one possible planetary system. One of the advantages of this 
technique is the direct measurement of the luminosity and 
projected orbits of the planets, making the physical charac- 
teri zation of the system and the individual planets possible. 

arois et ail l|2008h used a procedure for estimating the 
mass of the objects around HR 8799 in order to discriminate 
whether they were real planets or brown dwarfs (BD). This 
procedure can be applicable to any direct imaging detection, 
and it is based on the comparison of theoretical evolutionary 
tracks of BD and giant planets with observations in an Age - 
Luminosity diagram. Since luminosity is a direct observable, 
this technique is limited by the accuracy of the age determi- 
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nation and the theoretica l models used for this comparison 
l|Reidemeister et al.|[2009h . 

The A5V spectral type star HR8799 (V342 Peg, 
HP 218396, HIP 114189 ) has been extensively studied. 
ISchuster fc Niss cn (1986) firstly reported this star as a pos- 
sible SX Phoenicis object (these stars are puls ating subd- 
warf s tars with periods larger than one day). IZerbi et al.l 
(|l999h classified it as one of the 12 first 7 Doradus pulsators 
known. This pulsating stellar type is composed of Main Se- 
quence ( MS) stars in the low er part of the classical instabil- 
ity strip (|Dupret et al.|[2005T ). with periods around one day. 
That means that their pulsating modes are asymptotic g- 
modes a nd they are suitable to study the deep interior of 
the star. iGrav fc Kavel (|l999l ) obtained an optical spectrum 
of HR8799, and assigned an spectral type of kA5 hFO mA5 
V A Bootis (see that paper for the meaning of this specific 
nomenclature), reporting a metallicity of [M/H]= —0.47. 
The A Bootis nature of the star means that it has solar sur- 
face abundances of light elements, and subsolar abundances 
of heavy elements, the internal metallicity of the star being 
unknown. They also noted that HR8799 may be also a Vega- 



2 A. Moya et al. 



Table 1. Age of HR8799 in the literature. The methods used 
are: a) Stellar kinematics groups (proper motions), b) HR di- 
agram position compared with isochrones, c) HR diagram po- 
sition compared with models, d) IR excess. 



Age range 


Reference 


Method used 


in Myr 






[20, 150] 


Moor et al. (2006) 


a 


[30, 160] 


Marois et al. (20081 


a, b 


[50, 1128] 


Sone et al. (2001) 


c 


30 


Zuckerman & Sone (2004) 


a, b 


30 


Rhee et al. (2007) 


a, b 


[30, 730] 


Chen et al. (2006) 


a, b, d 



type star, characterized by a far IR excess due to a debris 
disk. Up to now, only three A Bootis stars h ave been reported 
to be 7 Doradus pulsators: HP 218427 jRodrfguez et al.l 
l2006al). HP 239276 JRodrfguez et al.ll2006bl ). and HR8799. 

iMarois et al.l l|2008l ) estimated the age of the star using 
four different methods (see Section 2 for details). The con- 
clusion of that work was that HR8799 is a young MS star 
with a n age in the range [30, 160] Myr. iReidemeister et al.l 
(2009) added another element for estimating the age of thi s 
star: the infrare d excess rat i o, stu died by ISu et al.l (|2006l ) 
and also used bv lChen et al.l l|2006l ). None of these determi- 
nations are conclusive, as we explain in Section 2. 

In this work, a comprehensive modelling of HR 8799 has 
been done in order to estimate its age and mass. The A Boo- 
tis nature of this star, and its 7 Poradus pulsations are the 
bases of the determinations presented here. All t he technical 
detai ls are presented in an accompanying paper (|Mova et al.l 
|2010| . Paper I), where, a complementary study on the A Boo- 
tis nature of this star is developed. 



2 PREVIOUS AGE DETERMINATIONS OF 
HR8799 

Most of the procedures to determine the age of isolated stars 
in the literature, such as activity, surface lithium abundance, 
or the position in the HR diagram, are useful for young 
and for low-mass stars, however these procedures cannot be 
acc urately applied to an A5 MS star. 

iKave and Strassmeierl l|l998l ) studied the chromospheric 
activity of some 7 Doradus stars, in particular HR 8799, us- 
ing the Ca II H&K lines, and found that all the stars studied 
have very low activity, HR8799 being the least active of all. 
This conclusion has an impact in its age determination since 
activity de creases with age in the early evolution ary stages 
of the star (Soderbl om. Duncan, fc John son 1991). The very 
low activity of HR 8799 suggests that it is a Main Sequence 
star. 

HR 8799 has a debris disk that has b een u s ed for an es- 
timation of its age. I Zuckerman fc Sond (|2004l ): Im oor et al] 
( 2006) studied the age of a number of stars with debris disks, 
searching for a relationship between the fractional IR lumi- 
nosity LiR/Lboi and the age. HR8799 was present in the 
samp les covered in these works ; they estimated a n age of 30 
Myr ijZuckerman fc Sondllobll: iRhee et al.ll2007l ) or in the 



range [20, 150] Myr ijMoor et al 



120061 ). Note that this age 



estimation was purely statistical, and has to be taken it with 



caution. In addition, the debris disks have cataclismic ori- 
gins (i.e. highly chaotic and unpredictable), and hardly can 
be r elated with the age of the star. 

Maroi s et al . (2008) used four methods for estimating 
the age of HR8799. The first one was the measurement of 
its radial velocity and proper motion (UVW). Comparing 
with other stars and associations with similar values of these 
quantities, they found an age range similar to that reported 
in previous works using also proper motions ([30, 160] Myr). 
They warned that this method is not always reliable. 

The second method was the position of the star in the 
HR diagram (for instance, see lPont fc E ver 2004). Compar- 
ing this position with isochrones of known stellar clusters 
one can infer upper limits to the age of the star. This ar- 
gument was used to confirm the age obtained with the sta- 
tistical UVW method. Nevertheless, the A Bootis nature of 
this star, and its unknown internal metalicity, makes the HR 
diagram position of t he star an inaccu rate method for the 
estimation of the star. ISong et al] l |200ll ) studied the age of a 
set of A-type stars (in particular HR 8799) using Stromgren 
photometry and isochrones, taking into account the stellar 
rotation. They concluded that, regarding the position in the 
HR diagram, the age estimated for this star is in the range 
[50,1128] Myr. 

The third method was based on the fact that A Boo- 
tis stars are generaly young. This argument is not correct. 
IPaunzen et al] ( 2002) found that most of the known A Bootis 
stars are stars between the ZAMS (zero-age main sequence) 
and the TAMS (terminal-age main sequence), with a mean 
age of 1 Gyr. 

The fourth method was based on the assumption that 
7 Doradus pulsators are probably young stars. There is no 
evidence for this assertion. 7 Doradus pulsations are origi- 
nated from the blocking of the radiative flux at the base of 
the outer convective zone (|Guzik e t al. 2000). Therefore, it 
is mainly a thermodinamic effect, dominated by the temper- 
ature of the star and the efficiency of the convective trans- 
port. The existence, or not, of 7 Doradus pulsations does 
not provide information about the age of the star, at least 
for a range in ages within the main sequence. 

Table 1 shows a summary of the different age determi- 
nations of this star found in the literature. 



3 ASTEROSEISMOLOGICAL ESTIMATION 
OF THE AGE OF HR8799 

In the previous section we indicated that there is no method 
that can clearly provide the age of the HR8799 planetary 
system. Of all the methods mentioned in the previous sec- 
tion, three seems to support a young age for our star. Of 
those, we will only consider those based on proper motions 
and on the position in the HR diagram, the third being an 
statistical method (viz., the luminosity of debris disks) with 
little accuracy in age determination. Our first step is to re- 
vise the physical parameters of the star and its position in 
the HR diagram with the aim of checking whether the over- 
lap in ages provided by this method and the proper motions 
is s upported 

iGrav fc Kavd (| 19991 ) obtained spectroscopic observa- 
tions of this star, providing accurate values for the stellar 
luminosity, T e g, and log g (see Table 2). The use of different 
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Table 2. Physical characteristics of HR8799 



Models of HR8799, no scismolottical constraints used 



T cS (K) 


7430±75 






logg (cms -2 ) 


4.35±0.05 


Grav & Kave 


(1999) 


L(L ) 


4.92±0.41 






vsini (km s —1 ) 


37.5±2 


Kayc and Strassmcicr 


(1998) 


it (mas) 


25.38±0.85 


van Lccuwcn 


(2007) 



bolometric corrections changes significantly the value of the 
luminosity of the star. We have u sed the Virtual observa- 
tory tool VOSA l|Bavo et al.ll2008l ) to avoide the necessity 
for estimating the stellar luminosity using bolometric cor- 
rections. Models with realistic metallicities best fitting the 
observations provide a luminosity only 0.1 L/Lq larger than 
that given in lGrav fc Kavd (|l999h and, therefore, within the 
errors. On the other hand, the use of different parallaxes in 
the literature changes the value of the absolute magnitude 
of this star by less than 0.1 L/Lq (see Paper I for details). 

We have developed a grid of equilibrium models ob - 
tained with the CESAM code l|Morel fc Lebretonl [2003) . 
with physics appropiate for main sequence A type stars (see 
Paper I) . In particular, the abun dance mixture used is that 
given in lGrevesse fc Noelsl (| 19931 ). To obtain the initial hy- 
drogen and helium abundances we have used a primordial 
helium content Yp = 0.235 and a mean enrichment law 
AY/AZ = 2.2. The use of other determinations o f these 
parameter in the literature (|Claret fc W illcms 20021. for ex- 
ample) does not change significantly the initial abundances 
used. The main approximation taken in the models with a 
possible influence in this study is the absence of updated 
internal chemical transport mechanisms in the equilibrium 
models. We vary the mass (in the range [1.25, 2.10] Mq with 
steps of 0.01 Mq), the metallicity (with values [M/H]=0.08, 
—0.12, —0.32, and —0.52), the Mixing-Length parameter 
MLT (values 0.5, 1, and 1.5), and the overshooting (values 
0.1, 0.2, and 0.3). The internal metallicity has been regarded 
as a free parameter due to the A Bootis nature of the star, 
which does not represent its inter nal abundances . The mass, 
estimated to be 1.47 ±0.3 M bv lGrav fc Kavd (|l999l) . has 
been also regarded as a free parameter since it has not been 
directly determined. 

If we search for models fulfilling the temperature, log g, 
luminosity and radius observed for the star, only a small 
amount of them remain, all in an age range [10, 2337] Myr 
(Fig. 1). A 18.1% of the models has the age estimated in 
the literature for this star ([30, 160] Myr). The gaps in the 
figure are due to the mesh of our grid of models. This result 
shows that the HR diagram position alone is not an accu- 
rate method to estimate the age of HR8799. The reason is 
the unknown internal metallicity of the star due to its A 
Bootis nature. Therefore, we need additional information to 
estimate the age of this planetary system. 

A comprehensive asteroseismic modelling of HR8799 
have been also developed. We use the presence of 7 Do- 
radus pulsations found in this star. As 7 Doradus stars are 
asymptotic g-mode pulsators, they are very good candidates 
for testing the internal structure of the star, in particular its 
internal metallicity and age, both having a large influence in 
the bouyancy rest oring force, measured through the Brunt- 
Vaisala frequency (|Mova et al . 2005; S mevers fc Movall2007l ; 
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Figure 1. T e g - Age diagram for the models of our grid fulfilling 
the spectroscopic observations, i.e., T^yy, logg and luminosity 
(see Table 2). The age range obtained is [10, 2337] Myr. 



iMiglio et alj liooil ). To do so, some of the most updated 
tools adapted for this purpose are used: 1) the evolution- 
ary code CESAM dMorel fc Lebretonll2008t). and 2) the pul- 
satio n code GraCo ( Mova et al.l |2004| ; iMova and Garridq 
120081) . T he latter has been tested in the work of the ESTA 
activity (|Mova et al.ll2008l ; iLebreton et al.l |200S| . Evolution 
and Asteroseismic Tools Activities). Using these tools, we 
performed a theoretical study of HR8799 in an attempt 
at constraining physical and theoretical parameters. In this 
work we will follow the same scheme used for the study of 
the 8 Scuti pulsators RV Arietis, 29 Cygn is and 9 Aurigae 
|Casas et al)l2006l . 120091 ; IMova et al.ll2006l ). 

A description of the tools and procedures used can be 
found in Paper I. In this work we focus only on the determi- 
nation of the mass and age of the star. These two quantities 
are the physical properties of the star with a larger impact 
in the understanding of the nature of the planetary system. 

An impo rtant unknown for the present study is the rota- 
tion velocity. iKave and Strassmeier] (1 19981 ) gave a projected 
rotation velocity for this star of vsini = 37.5 ± 2.0kms _1 . 
Therefore, an accurate determination of v moves to an ac- 
curate determination of the visual angle i. Several works 
in the literature tried to provide estimations for i with no 
conclusive results. This is an important condition for the 
modelling. 

In the present study, we use the F reque ncy Ratio 
Method (FRM) developed in IMova et all (|2005l ) (see Sec- 
tion 3.3 of Paper I) for the asteroseismological modelling 
of this star. The FRM is based on the analytical descrip- 
tion of the frequencies in the asymptotic regime. This 
procedure has a limited range of validity in stellar r ota- 
tion velocities, as demonstrated in ISuarez et al.l l|2005h . In 
that paper, the authors found a maximum rotation veloc- 
ity for the correct application of the FRM, for a standard 
7 Doradus star, of around 60k ms~ . On the other hand, 
iTurcotte fc Charbonneaul (|l993h found that meridional cir- 
culation cannot destroy the accretion pattern (in the ac- 
cretion/diffusion scenario, based on the accretion of inter- 
stellar medium by the star, and the mixing of these elements 
in the stellar surface due to diffusion and rotationally mix- 
ing processes) for a rotation velocity below 125kms _1 . For 
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Table 3. Acceptable models depending on the physical constraints. In the less favorable case, the first (second) mass range 
is linked with the first (second) age range shown. 



Constraint 


Mass 


Age 


% of models with [30, 160] 




in Mq 


in Myr 


in Myr 


HR position 


[1.25, 1.27], [1.32,1.35], [1.40,1.48] 


[10, 2337] 


18.1 


Complete procedure 


1.32 


[1126, 1486] 


0.0 


Comp. proc. (less favorable case) 


[1.32, 1.33], [1.44, 1.45] 


[1123, 1625], [26,430] 


16.7 



a larger rotation velocity, theory cannot ensure that the ac- 
creted abundances can remain in the stellar surface enough 
time to be observed. These values provide two lower limits 
for the inclination angle. 

Obtaining the rotation velocity as a function of i for any 
stellar radius in the observed range we found that i > 18° 
would be the requirement imposed by the A Bootis nature of 
HR8799, and i > 36° would be the requirement for the FRM 
to be applicable. We want to point out that if the inclination 
angle is between both values, our analysis would be possible 
but inaccurate, since we cannot ensure that the real solution 
is part of that provided by the FRM (|Suarez et al.ll2005r i. 

Our first study was done for a moderate rotation case 
in the Frequency Ratio Method. This means an inclination 
angle i w 50° and a rotation velocity V I0 t ~ 45kms _1 . In 
this case, the complete procedure is very discriminant and 
only models with M= 1.32 Mq fulfill all the observational 
constraints. The acceptable age range is [1126, 1486] Myr, 
far fro m the value estimated for this star by iMarois et al.l 
(|2008t ) ([30,160] Myr). In Fig. 2 (left panel), the acceptable 
set of models for this rotation velocity is shown in a T c ft-Age 
diagram. The T e g axis limits are provided by the spectro- 
scopic observations. 

We have also studied the less favourable case in the 
FRM, i.e. an inclination angle i « 36°, that is, V TO t ~ 
60kms _1 . In this case, the complete procedure is not com- 
pletely discriminant and new models fulfill all the obser- 
vational constraints. Two mass ranges are now acceptable 
(M=[1.32, 1.33] and [1.44, 1.45] M Q ), and the compatible 
age range is divided into two blocks: [1123, 1625] Myr and 
[26, 430] Myr respectively (Fig. 2, right panel). The two 
models with ages around 800 Myr have been accepted due 
to the conservative application of the instability analysis (see 
Paper I), but it their belonging to the set of acceptable mod- 
els in unlikely. In this less favourable case, the percentage of 
accepted models with ages in the range commonly used for 
this star is 16.7%. A summary of the results presented here 
is found in Table 3. 



4 CONCLUSIONS AND FUTURE 
PERSPECTIVES 

In the present work, an analysis of the age determination of 
the planetary system HR8799 has been done. The results 
found in the literature are not conclusive, and the only valid 
argument to estimate the age of the star is that using its 
radial velocity and proper motion, but it is an estatistical 
argument needed of additional estimations. 

The only valid argument used to estimate the age of 
the star is its radial velocity and proper motion, but it is a 
statistical argument needed of additional estimations. The 



main complementary argument is the position of the star in 
the HR diagram. In this work we have demostrated that this 
procedure does not provide accurate age estimations due to 
the A Bootis nature of HR8799. This nature hides the real 
internal metallicity of the star. The main consequence of this 
result is that the models fulfilling observations are in a range 
of ages [10,2337] Myr, a much broader range that one est i- 
mated by other authors of [30,160] Myr (|Marois et al.ll2008h . 
Only a small amount (18.1%) of models in our representative 
grid have ages in the range claimed in the literature. 

Therefore we need aditional constraints for an accurate 
estimation of the age and mass of the star (these two quan- 
tities are the physical characteristics of the star with larger 
impact in the understanding of the planetary system). We 
have taken advantage of the 7 Doradus pulsations of the star 
to better estimate these values with the help of asteroseis- 
mology. 

A comprehensive asteroseismological study of this star 
has been develop ed. This study is described in detail in 
(|Mova et alj|20ld , Paper I). The main source of uncertainty 
of the procedure is the unknown rotation velocity of the star. 
We have analysed the possible results depending on the in- 
clination angle i. There is a range of angles where this study 
is not accurate (i — [18°, 36°]). 

For angles around i = 36°, the models fulfilling all the 
observational constraints have masses in two separate ranges 
of M=[1.32, 1.33], [1.44, 1.45] Mq. The age of the system 
is constrained in two separate ranges: [1123, 1625] Myr and 
[26, 430] Myr respectively. A percentage of 16.7% of the 
models are in the range given in the literature, i.e., young 
ages. A consequence of this result is that, in the case of the 
youngest age range, the predicted masses of the observed 
planets are [5,14] Mj up for the most luminous planets, and 
[3,13] Mj up for the less luminous one. The oldest age range 
predicts masses for t he three objects in the brown dwarfs 
domain (see Fig. 4 of lMarois et al.ll200sl ). 

In the most favourable case for the procedure used in 
this work, i.e. inclination angles of around i ~ 50°, astero- 
seismology is very accurate, and the star would have an age 
in the range [1126, 1486] Myr, and a mass M=1.32 Mq. This 
age range implies that the observed objects orb iting HR 8799 
woul d be brown dwarfs (following Fig. 4 in M arois et al.l 
l200St) . 

The lack of an accurate determination of the inclination 
angle is the main source of uncertainty of the present study. 
This angle has not been unambiguously obtained up to now, 
and its value would say whether the results of this study are 
actually applicable, and then, it can provide a very accurate 
determination of the age and mass of the star. In any case, 
one of the main conclusions of our study is that the range of 
ages assigned to this star in the literature is unlikely to be 
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Figure 2. T e g - Age diagram for the models of our grid fulfilling the spectroscopic observations plus all the asteroseismological constraints 
(FRM + multicolour photometry + instability analysis, see Paper I) for a rotation velocity Vrot ~ 45km s -1 (left panel) and V ro t ~ 
eOkrns" 1 (right panel). The age ranges obtained are [1123, 1625] Myr and [1123, 1625], [26,430] Myr, respectively 



the correct one. Only a stellar luminosity larger than that 
reported would allow young models with solar metallicity to 
fulfill all the observational constraints. 
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